Intermodal entanglement in Raman processes 
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The operator solution of a completely quantum mechanical Hamiltonian of the Raman processes 
is used here to investigate the possibility of obtaining intermodal entanglement between different 
modes involved in the Raman processes (e.g. pump mode, Stokes mode, vibration (phonon) mode 
and anti-Stokes mode). Intermodal entanglement is reported between a) pump mode and anti-Stokes 
mode, b) pump mode and vibration (phonon) mode c) Stokes mode and vibration phonon mode, 
d) Stokes mode and anti-stokes mode in the stimulated Raman processes for the variation of the 
phase angle of complex eigenvalue ai of pump mode a. Some incidents of intermodal entanglement 
in the spontaneous and the partially spontaneous Raman processes are also reported. Further it 
is shown that the specific choice of coupling constants may produce genuine entanglement among 
Stokes mode, anti-Stokes mode and vibration-phonon mode. It is also shown that the two mode 
entanglement not identified by Duan's criterion may be identified by Hillery-Zubairy criteria. It 
is further shown that intermodal entanglement, intermodal antibunching and intermodal squeezing 
are independent phenomena. 



I. INTRODUCTION 

Entanglement is one of the most important resources 
for quantum communication and quantum information 
processing. For example, it is well known that entangle- 
ment is essential for teleportation, dense coding, quan- 
tum information splitting etc. Thus we need entan- 
gled states to perform various important tasks related 
to quantum information theory. To do so, first we need 
a protocol to check, whether a state generally mixed is 
entangled or not? This is a very important issue in quan- 
tum information science and several inseparability cri- 
teria have been proposed for this purpose and ref- 
erences therein). In 1996, Peres [2| proposed the first 
inseparability criterion based on negative eigenvalues of 
partial transpose of the composite density operator. This 
criterion is sufficient and necessary for the detection of 
entanglement in (2x2) and (2x3) dimensional states, but 
is not necessary for higher dimensional states (see Q 
and references therein). Since the pioneering work of 
Peres, several other inseparability inequalities have been 
reported for two mode and multi-mode states 
Most of these criteria only provide sufficient condition 
of inseparability. Further, these criteria may be classi- 
fied into two sets [loj : A) set of criteria which cannot be 
directly tested through experiments and B) set of 

criteria which can be tested experimentally 
Experimentally testable inequalities involve variance or 
higher order moments of some observables. Since the ex- 
pectation values of physical observables can be measured 
experimentally, these set of inseparability criteria can be 
tested experimentally. 

The aim of the present work is not to study the in- 



separability criteria in detail but to study the possibil- 
ity of generation of multi-partite entangled state in two- 
photon stimulated Raman processes, as depicted in Fig. 
[1] The scheme is essentially a sequential double Raman 
process that can produce Stokes and antiStokes photons 
that show highly nonclassical correlation 11 ll and macro- 
scopic entanglement in two-photon laser |l2l |. To study 
the two-mode entanglement in the two-photon Raman 
processes it would be reasonable to use three criteria from 
set B. To be precise, we have chosen the two criteria of 
Hillery and Zubairy Q, Q and the criterion of Duan 
et al. [3j. Since all these three criteria are only suffi- 
cient, a particular criterion can detect only a subset of 
all sets of entangled states. Consequently, application 
of a single criterion may yield incomplete result. This 
is why we have used three experimentally testable in- 
separability criteria for our investigation of intermodal 
entanglement in stimulated, spontaneous and partially 
spontaneous Raman processes. 

Nonclassical properties of these Raman processes have 
been extensively studied. Initial studies were restricted 
to the short-time approximation [3|-[lB|. But recently 
some of the present authors have reported different non- 
classical effects (such as squeezing, antibunching, inter- 
modal antibunching and sub-shot noise photon number 
correlation) in stimulated and spontaneous Raman pro- 
cesses [13]- [13 without using traditional short-time ap- 
proximation technique. Our solution of Raman pro- 
cesses, which does not involve short-time approxima- 
tion, is found to reveal many facets of nonclassical ef- 
fects which were undetected by short-time approxima- 
tion technique. However, the possibility of observing in- 
termodal entanglement is not rigorously studied so far. 
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This fact has motivated us to study the intermodal en- 
tanglement in the double Raman processes. The present 
investigation is relevant for quantum communication for 
two reasons: Firstly because entanglement is an essen- 
tial resource for quantum communication and secondly 
because spontaneous Raman process is reported to be 
useful in the realization of quantum repeaters [2ll]-[22| 
which has its application in long distance high-fidelity 
quantum communication. 
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ical treatment where all four modes are considered quan- 
tum mechanical. If we look closely into the methodology 
adopted in the earlier studies we would quickly find that 
the approach adopted in the present paper is simpler and 
easily extensible to the other physical systems which are 
described by bosonic Hamiltonians. 

The paper is organized as follows. In the Section [XT] 
we have described the Hamiltonian of spontaneous and 
stimulated Raman processes and its operator solution. 
In Section Mil we have used the solution to show that 
it is possible to observe intermodal entanglement in Ra- 
man processes. The inseparability criteria used for this 
purpose are also described in this section. Finally, Sec- 
tion IIVI is dedicated to conclusions where we have briefly 
summarized the result of the present study and have also 
discussed the mutual relations among different nonclas- 
sical phenomena observed in the Raman processes. 



II. MODEL HAMILTONIAN 



Figure 1: (Color online) Two-photon stimulated Raman 
scheme. The pump photon is converted into a Stokes 
photon and a phonon. The pump photon can also mix 
with a phonon to produce an anti-Stokes photon. 

Here it is worthy to note that V Pefinova et al. (23| 
have recently studied the possibility of observing entan- 
glement in Raman process using the method of invariant 
subspace. They have followed an independent approach 
and have numerically computed the time dependence of 
a measure of entanglement. Earlier S V Kuznetsov et 
al. (24| studied the entanglement in the stimulated Ra- 
man process considering only two modes (Stokes mode 
and phonon mode) and taking the pump mode as the 
classical light source. Naturally Kuznetsov et al.'s work 
illustrated an incomplete scenario and failed to observe 
intermodal entanglement involving anti-Stokes mode and 
pump mode. To circumvent this limitation we have 
used here a completely quantum mechanical Hamilto- 
nian. Further, Pathak, Kfepelka and Pefina [25[ have 
recently investigated the possibilities of observing inter- 
modal entanglement in the Raman processes using the 
same Hamiltonian but with a short-time approximated 
solution. Their work is restricted by the intrinsic lim- 
itations of the short-time approximation. Such limita- 
tions may be circumvented by the analytical methods 
developed by us in recent past to study the stimulated 
Raman scheme [l7]]-[2(|. Those methods are systemati- 
cally used here and a relatively complete scenario of in- 
termodal entanglement in Raman processes is presented 
here. Interestingly, we have observed intermodal entan- 
glement between i) pump mode and anti-Stokes mode, 
and ii) Stokes mode and anti-stokes mode. These two 
intermodal entanglement was not observed in the earlier 
analytic studies [2J, l25fl . The beauty of the present study 
lies in the fact that analytic expressions for separability 
criterion are obtained by a completely quantum mechan- 



The Hamiltonian [lj|-[21|, |26j of our interest is 

H = uj a <r a + LO^b'b + uj c c'c + iv<i<v d 
+ g (atfrf + h.c.) + x (aatf + h.c.) , 



(1) 



where h.c. stands for the Hermitian conjugate. Through- 
out the present paper, we use h = 1. The annihila- 
tion (creation) operators a(aJ), b(w), c (cM , d(cft) corre- 
spond to the laser (pump) mode, Stokes mode, vibration 
(phonon) mode and anti-Stokes mode, respectively. They 
obey the well-known boson commutation relations. The 
quantities u> a , wt,, uj c and ujd correspond to the frequen- 
cies of pump mode a, Stokes mode 6, vibration (phonon) 
mode c and anti-Stokes mode d, respectively. The pa- 
rameters g and \ are the Stokes and anti-Stokes coupling 
constants, respectively. Coupling constant g (x) denotes 
the strength of coupling between the Stokes (anti-Stokes) 
mode and the vibrational (phonon) mode and depends 
on the actual interaction mechanism. The dimension of 
g and x are that of frequency and consequently gt and 
\t are dimensionless. Further, gt and %t are very small 
compared to unity. In order to study the possibility of in- 
termodal entanglement, we need simultaneous solutions 
of the following Heisenberg operator equations of motion 
for various field operators 



a = — i (uj a a + gbc + x°d) 

b = —i (ojf,b + gad) 

c = —i (cj c c + gab^ + x a d) 

d = — (cudd + x ac ) ■ 



(2) 



The above set of coupled nonlinear differential operator 
equations © are not exactly solvable in the closed ana- 
lytical form under weak pump condition. But when the 
pump is very strong the operator a can be replaced by 
a c-number and the above set of equations ^ can be 
solved exactly [16]. In order to solve these equations 
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under weak pump approximation we use the perturba- 
tive approach. Our solutions are more general than the 
well-known short-time approximation. Details of the cal- 



culations are given in our previous papers Here 
we just note that under weak pump approximation, the 
solutions of Eq. ([2} assume the following form: 



a(t) = /ia(0) + / 2 6(0)c(0) + / 3 ct(0)d(0) + / 4 a t (0)6(0)d(0) + f 5 a(0)b(0)tf (0) 

+ / 6 a(0)ct(0)c(0) + / 7 a(0)ct(0)c(0) + / 8 o(0)d+(0)d(0) 

b(t) = 3l fo(0) + . 92 a(0)ct(0) +53^(0)^(0) + 54 c t2 (0)d(0) + . 95 &(0)c(0) C t(0) 

+ 5r 6 &(0)a(0)at(0) 

c(t) = h lC {0) + /i 2 a(0)6t(0) + ft 3 ot(0)d(0) + M f (0)ct(0)d(0) + /i 5 c(0)a(0)ot(0) ' 

+ &8c(0)&(0)&t(0) + /i 7 c(0)d+(0)d(0) + /i 8 c(0)at(0)a(0) 

d(f) = iid(0) + Z 2 a(0)c(0) + / 3 a 2 (0)&t(0) + Z 4 6(0)c 2 (0) + Z 5 ct(0)c(0)d(0) 

+ leo(0)o+(0)d(0) 



The functions /,;, <?j, ft,,; and Z; are evaluated from the 
dynamics under the initial conditions. In order to apply 
the boundary condition, we put t = 0, in the first term 
of the Eq. ©. It is clear that /i(0) = gi (0) = hi(0) = 
h(0) = 1 and MO) = 9i (0) = hi(0) = h(0) = (for 
i = 2, 3, 4, 5, 6, 7 and 8). Under these initial conditions 
the corresponding solutions for /j(t), <7i(t), hi(t) and Zj(i) 
are obtained as given in the Appendix. 

The solutions Eqs. ©, ([XT ]) -([X3 ]) are valid up to the 
second orders in g and x- Interestingly, there is no re- 
striction on time t. For example, / 2 rises indefinitely with 
the increase of time t. Clearly, the divergent nature of the 
parameters /j, gi, hi and Zj become more explicit as the 
order of the perturbation theory is increased. The sec- 



ular nature is a direct outcome of the perturbation the- 
ory. In the present investigation the secular term is not a 
problem since we consider small interaction time. Small 
interaction time also ensures that the damping term con- 
tributes insignificantly. Here Awi = w& + uj c — u a and 
Aw 2 = Lj a + uj c — LJd ■ Normally, the detunings Awi and 
Aw 2 are extremely small. In the present investigation, 
we however assume that the small (non-zero) detuning 
is present and hence Awi ^ and Aw 2 ^ 0. Here we 
have used |Acji| = 0.1 MHz and |Aw 2 | = 0.19 MHz. Of 
course, in Eqs. (| Al|) - (| A4|) we have neglected the terms 
beyond the second order in g and \- Now we may use 
Eq. ([5]) to obtain the temporal evolution of the number 
operators of various modes as 



N a (t) = |M 2 at(0)a(0) + |/ 2 | 2 &t(0) C t(0)&(0)c(0) + |/ 3 | 2 c(0)d(0)t c t(0)d(0) 

+ [/r/2a f (0)&(0)c(0)+ /*/ 3 at(0)ct(0)d(0) + / 1 */ 4 a t (0)a t (0)6(0)d(0) 

+ f?fc(aH0)a(0) + at(0)a(0)&t(0)6(0)) + / 1 V 6 a t (0)a(0)ct(0) C (0) 

+ / 1 */ 7 a t (0)a(0)ct(0)c(0) + A*/ 8 at(0)a(0)dt(0)d(0) 

+ /|/ 3 & t (0)ct 2 (0)d(0)+h.c.l , 



N b (t) = |. gi | 2 6t(0)6(0) + | 32 | 2 at(0)c(0)a(0)ct(0)+ [ 5l32 6t(0) a (0)ct(0) 

+ ffl * ff3 6t(0)a 2 (0)dt(0) + 5l * 54 6t(0) c t 2 (0)d(0) + (^5 (6* (0)6(0) 

+ 6t(o)6(o) c t(0)c(0)) + <7iffe(& + (0)6(0) + &t(0)6(0)at(0)a(0)) 

+ h.c], 



N c (t) = |/n| 2 ct(0)c(0) + |ft 2 | 2 at(0)6(0)a(0)6t(0) + |/i 3 | 2 a(0)dt(0)at(0)d(0) 

+ [/ij/i 2 cf (0)o(0)6+(0) + hlh 3 c^0)at(p)d(0) 

+ /it/i 4 ct(0)&+(0)c+(0)d(0) + hlh 5 (cH0)c(0) + ct(0)c(0)o+ (0)o(0)) 

+ fe;ha( c t(0)c(0) + c+(0)c(0)&t (0)6(0)) + ft^ 7 c t (0)c(0)dt(0)d(0) 

+ ^ft 8 ct(0)c(0)at(0)a(0) + /i^ 3 at 2 (0)6(0)d(0)+h.c.l , 



and 
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N d (t) = |/i| 2 dt(0)d(0) + |/ 2 | 2 at(0)a(0)ct(0)c(0)+ [/^ 2 a(0)c(0)dt(0) 

+ Z^3a 2 (0)6t(0)dt(0) + Z^4&(0)c 2 (0)dt(0) + ;*; 5 ct(0)c(0)dt(0)d(0) . (7) 
+ ZJi 6 dt(0)d(0)(l + at(0)o(0)) + h.c] . 

I 



In the following section these number operators will be 
used to study the intermodal entanglement in Raman 
processes. 

III. INTERMODAL ENTANGLEMENT 

In order to investigate the intermodal entanglement for 
various coupled modes, we assume that all photon and 
phonon modes are initially coherent. In other words, the 
composite boson field consisting of photons and phonon 
are in the initial coherent state. Therefore, the composite 
coherent state arises from the product of the coherent 
states \ai), {012)1 l a 3): ail d \ola) which are the eigenkets 
of a, b, c and d respectively. Thus the initial composite 
state is 

1^(0)) = <g> \a 2 ) ® \a 3 ) ® |o4>. (8) 

It is clear that the initial state is separable. Now the field 
operator a{t) operating on such a multi-mode coherent 
state gives rise to the complex eigenvalue ot\{t). Hence 
we have, 

a(0)|V(0)} = ai|V(0)>, (9) 

where \a\\ 2 is the number of input photons in the pump 
mode a. In a similar fashion we have three more complex 
amplitudes ct2(t), as(t) and 014, (t) corresponding to the 
Stokes, vibrational (phonon) and anti-Stokes field mode 
operators b, c and d respectively. Clearly, for a sponta- 
neous process, the complex amplitudes are a 2 = «3 = 
Q'4 = and ai 7^ 0. For partial spontaneous process, 
the complex amplitude a\ and any one of the remaining 
three eigenvalues are not equal to zero while the other 
two complex amplitudes are zero. On the other hand, 
for a stimulated process, the complex amplitudes are not 
necessarily zero. In our present investigation we consider 
ai = \ai \ er % ^ and the other eigenvalues for the Stokes, 
vibrational (phonon) and anti-Stokes field modes are real. 
The aim of the present work is to investigate the pos- 
sibility of intermodal entanglement in the spontaneous, 
partially spontaneous and stimulated Raman processes. 
To do so let us begin with the investigation of two mode 
entanglement using Hillery and Zubairy's criteria. 

A. Two mode entanglement 

There are two criteria due to Hillery and Zubairy 0|- 
@| . The first one is 

(N a N b ) - \(atf)\ 2 < 0. (10) 



On the other hand, the second criterion is given by 

(N a )(N b )-\(ab)\ 2 <0. (11) 

From here onward we will refer to these criteria as HZ- 
1 and HZ-2 criterion respectively. In addition to these 
two criteria, we will also use the Duan's inseparability 
criterion due to Duan et al. Q : 

(Aa t Aa)(A6 t A6) - (AaA6)| 2 < 0. (12) 

In the criteria Eqs. (fT0)) - (fT2|) . a and b are annihilation 
operators for two arbitrary modes. They are not limited 
to the pump mode and the Stokes mode only. 

We note that all the above criteria are only sufficient 
(not necessary) for detection of entanglement. Keeping 
this fact in mind, we have applied all these criteria to 
study intermodal entanglement between different modes 
of Raman Hamiltonian and have observed intermodal en- 
tanglement in various situations. 

Let us first investigate the possibility of two mode en- 
tanglement in Raman process using HZ-1 criterion. From 
Eqs. ©, 0, © and © we obtain 

(N a N b )- |(afet)| 2 = |/ 3 | 2 |« 2 | 2 |a 4 | 2 

+ |.92| 2 (jai| 4 ~ \oi\ 2 |a 2 | 2 ) • 

(13) 

Consequently, for spontaneous Raman process Eq. (fT"3")) 
reduces to 

(N a N b )- |( a fet)| 2 = |g 2 | 2 | ai | 4 (14) 

It is evidently clear that the right hand side of the 
Eq. (|14p is always positive. Hence HZ-1 criterion does 
not show any signature of intermodal entanglement be- 
tween pump and Stokes modes in the spontaneous Ra- 
man process. For partially spontaneous Raman process 
(|oi| ^ 0, |ata I ^ 0, |a 3 | = |a 4 | = 0), the Eq. Q3) re- 
duces to 

(N a N b ) - |(afet)| 2 = | 52 | 2 (| Ql | 2 - |a 2 | 2 ) K| 2 . (15) 

It is clear that the entanglement is possible in the par- 
tially spontaneous Raman process only when |a 2 | > 
\a\\ , i.e. the number of Stokes photon is more than the 
number of pump photons, which is not the usual case. 
According to the HZ-1 criterion of Eq. (fT0|) , it is clear 
that the negative values on the right hand side of Eq. 
(fT"3")) would indicate the presence of intermodal entangle- 
ment between the pump mode and the Stokes mode in 
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stimulated Raman process. To investigate the possibil- 
ity of intermodal entanglement in the stimulated Raman 
process we have used \ = 9 = 10 5 Hz, |ai| = 10, |ck2 | = 8, 
|a 3 | =0.01 and [au| = 1 [27j|. We have plotted the right 
hand side of ([TBI in Fig. [2^ which does not show any sig- 
nature of intermodal entanglement between pump mode 
and the Stokes mode in the stimulated Raman process. 

Here we would like to note that once we have an an- 
alytic expression for the HZ-1 or HZ-2 or Duan crite- 



ria in stimulated Raman process, it is straightforward to 
study the special cases of: i) spontaneous Raman pro- 
cess, where a 2 = a 3 = a 4 = but a\ ^ and ii) partial 
spontaneous Raman process, where ot\ ^ and any one 
of the other three a» (i = 2, 3, 4) is non-zero. 

The same technique used in the above case is now 
adopted to obtain the following equations for the study 
of intermodal entanglement in stimulated Raman process 
using HZ-1 criterion: 



(N b N c ) - | (6c t ) | 2 = \g 2 \ 2 (z |«i | 2 |a 3 | 2 + 3 |ai| 2 \a 2 \ 2 + |«i| 2 - \a 2 \ 2 |a 3 | 2 ) + \h 3 \ 2 \a 2 \ 2 |a 4 | 2 ^ 
+ [{hXh^ict^at, + 2glgia 2 a\a% + h 2 hlala 2 al] + cc] 

(N a N d ) - \(aS)\ 2 = |/ 3 | 2 (|a 3 | 2 + M 4 + l«i| 2 K| 2 - K| 2 |a 4 | 2 ) - \l 2 \ 2 (|a 3 | 2 + K| 2 |a 3 | 2 ) (17) 
(N b N d ) - |(Mt)| 2 = \g 2 \ 2 \ ai \ 2 \a 4 \ 2 + [{Itkaja^al + c.c.}] (18) 
(N c N d ) - |(crft)| 2 = \l 2 \ 2 ( 2 | ai | 2 + 2 M 2 |a 4 | 2 - 2 |a 4 | 2 - |« 4 | 4 - |a 3 | 2 |a 4 | 2 ) + \h 2 \ 2 \ ai \ 2 |a 4 | 2 (19) 

(AUV C ) - |(act)| 2 = |/ 2 | 2 ^2 K| 2 + |a 4 | 4 + M 2 |a 3 | 2 - 4 \a 2 \ 2 - 2 | ai | 2 |a 2 | 2 - 2 \a 2 \ 2 |« 3 | 2 ) 

+ l/sl 2 (l« 4 | 2 + 3 |a 3 | 2 M 2 + 3 |a 4 | 2 |a 4 | 2 - K| 2 |a 3 | 2 ) + [{ftf 3 ata* 3 a 4 ( 2 0) 
+ h^h^a* 2 : a 2 a 4 + ,f 2 .f3a 2 a 3 2 a 4 ) + cc] . 



The right hand sides (RHS) of Eqs. (fT6 ]) - (j2"0]l are plot- 
ted in Fig. - Fig. [2f. It is interesting to note that 
the presence of intermodal entanglement in stimulated 
Raman process is observed between i) the Stokes mode 
and the vibration (phonon) mode (Fig. [5})), ii) the pump 
mode and the anti-Stokes mode (Fig. [5t), hi) the Stokes 
mode and the anti-Stokes mode (Fig. [5Ji) and iv) the 
pump mode and vibration mode (Fig. [5f). However, no 
signature of intermodal entanglement is observed in the 
other two cases (Figs. [5^ and e). Further, it does not 
show the presence of genuine entanglement between any 



three modes of the system. Interestingly, with the suit- 
able choice of the complex eigenvalues a,; it is possible 
to observe the signature of the intermodal entanglement 
using HZ-1 criteria in partially spontaneous Raman pro- 
cess in several ways but no such signature is observed for 
the completely spontaneous Raman process. 

Since the HZ-1 criterion is only sufficient, we might 
have failed to detect some intermodal entanglement. In 
an attempt to detect such intermodal entanglement using 
HZ-2 criterion (TT]), we have used Eqs. ©, ©-jZl) and 
© to yield: 



(N a )(N b )-\(ab)\ 2 = 



and 



I52I 2 Kl 4 

/l/2Sl32) |«l| 2 |«2| 2 



I/3I \a 2 \ 2 \a 4 \ 2 - [(glge 



(21) 



(N b )(N c )-\(bc)\ 2 = 



|32| 2 Kl 2 |a 3 | 2 



[(h\h 



2 a\a 2 a 3 



-\h 2 f 
(hiht 



1 



«2 



|ft, 3 | 2 \a 2 \ 2 |a 4 | 2 



gig 2 hih 3 ) a 2 a 3 a\ + h 2 h 3 a\ a 2 a 4 



(22) 



K^s « 2 | 2 |a 3 | 2 +g 4 g 2 hlh 2 |ai| 2 |a 3 ' 2 



r +c.c. 



(i 



(N a )(N d )-\(adW 



l/ 3 | 2 KI 4 



6 |ai| 2 |a 4 | 2 + c.c 



(23) 



(N b ) (N d )-\(bd}\< 



M 2 |ai| 2 \a 4 \ 2 



3 af a 2 a! 4 + c.cl 



(24) 



(JV C ) (N d )-\(cd)\< 



\h 2 \ 2 lail 2 |a 4 | 2 



l^3| 2 |Q4| 4 



|ct 3 | |a 4 | 



(25) 



(N a )(N c )-\(acW 



|/i 2 | 2 |a<i| 4 



|/13 
2 



| 2 ( 



«4 



lad |a 4 | 2 



|/ 3 | 2 |a 3 | 2 |a 4 | 2 



[(/i*ft. 3 a*a3a 4 



/i*/i 8 |ai| |a 3 | + h^h^al 2 a 2 a 4 - h{h 5 \ai\ \a 3 
fif2hlh 1 a 2 alal + f*hh%h\ \a 3 \ 2 |a 4 | 2 ) + c.c. 



(26) 



r 



From the closed form analytic expressions Eqs. (|21l) - 
|), it is possible to obtain the signature of intermodal 
entanglement in various cases. Interestingly, we obtain 
intermodal entanglement between Stokes mode and vi- 
bration mode for spontaneous Raman process. The in- 
termodal entanglement Eqs. (f2"T j) - (|2l))) for stimulated Ra- 
man processes are illustrated in the Fig. [3£i-Fig. |3f. In 
accordance to HZ-2 criterion, negative values of the ordi- 
nates indicate the signature of entanglement. Therefore, 
the intermodal entanglement is observed in: i) Stokes 
mode and vibration mode, ii) Stokes mode and anti- 
Stokes mode and iii) pump mode and vibration mode. 
However, there is no signature of intermodal entangle- 
ment in the remaining cases for Stimulated Raman pro- 
cesses. It is possible to obtain the intermodal entangle- 
ment for various partially spontaneous Raman processes. 
However, these results are not exhibited in the present 
text. It is interesting to note that HZ-2 criterion failed to 
detect intermodal entanglement between, pump and anti- 
Stokes mode. Thus the Raman process provides a very 
nice example of physical system where it can be shown 
with physical example that these inseparability criteria 
are only sufficient. Still there are two situations where 
we have not found the signature of intermodal entangle- 
ment. Let us see what happens when we apply another 
sufficient but not necessary criterion of inseparability. 

Now Duan criterion Eq. (fT2f for the intermodal en- 
tanglement can also be written as [28| 



D = ((Aw)' 



(Av) 



-2 < 



where 



and 



l 

i\/2 



(6 + 6t)] 



[(o_at) + (6-6t)] 



(27) 



(28) 



(29) 



Using Eqs. ((3]), © and (J2T1) - (J29J) we can obtain analytic 
expression for D the left hand side of Duan et al. criterion 
Eq. (H7D, 

For mode a and b, 



D ab = 



2[|/ 3 ||04 



S2I 2 K| 2 



(30) 



(2/ifl3 + Ugl + hgt) a\a A + c.c.}] 



for mode a and c 



D nr = 2 



I/3I |ct 4 | 2 + \h 2 \ 2 |ai| 2 + \h 3 \ 2 |a 4 | 2 



+ M* + Ms) a i a a} + c - c -l > 
for mode b and c 

D bc = 2[\g 2 \ 2 \ ai \ 2 + \h 2 \ 2 \a 1 \ 2 + \h 3 \ 2 \a 4 \ 2 
+ \ {(9i h 6 + 9sh* + g 2 h* 2 ) a^a 2 } + c.c] 

for mode a and d 



(31) 



(32) 



D„ 



\h\\^\ 2 + U(h l *2+m)a* 3 



+ (2/iiS + Uq + hl* 2 ) a\a 2 + f 8 l* iai a* 4 + c.c.}] , 

(33) 

for mode c and d 



;i| 2 + \h 3 \ 2 |a 4 | 2 ■ 



D cd = 2 [\h 2 \ 2 \ ai 

+ ±{(2hhl + hhZ + hhl)a2a s ^ 
+ (I5I1I + M7) a 3 a 4 + +c.c.}] , 



and for mode b and d 
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nc 
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Table I: Relation between different nonclassical phenomena observed in stimulated Raman process. Here nc stands 

for non-conclusive. . 



Right hand sides of equations (f5U |) - (j55)) are plotted 
in Fig. Fig 0J. It is clear that the intermodal- 

entanglement is observed only between the pump mode 



and anti-Stokes mode. The criterion is non-conclusive in 
all other cases. This is so because Duan criterion is only 
sufficient. 



IV. CONCLUSIONS 

We have clearly established that the stimulated Raman 
process can produce intermodal entanglement. The ob- 
servations are summarized in the Table |H Here it would 
be apt to note that recently Pathak, Kfepelka and Pefina 
(25| have investigated the possibilities of observing in- 
termodal entanglement in the Raman processes using a 
approximated short-time solution. They have identified 
intermodal entanglement in pump-phonon ac and Stokes- 
phonon be modes only. However, in addition to those 
two modes we have observed intermodal entanglement 
in pump-antiStokes ad and Stokes- antiStokes bd modes 
too. In addition to these, we explore the various pos- 
sibilities of getting intermodal entanglement in partial 
spontaneous Raman processes too. In this way, our so- 
lution is found more powerful compared to those of the 
solutions of Raman processes under short-time approxi- 
mation. Further, the use of short-time solution led to the 
monotonic nature of entan glem ent parameter as seen in 
Eqs. (11) and (12) of ref. [25||. As our solution is valid 
for all times and hence the entanglement parameters are 
free from this particular problem which is generally a 
characteristic of short-time solutions. Another earlier ef- 
fort to study the intermodal entanglement in the Raman 
processes by S. V. Kuznetsov [24| was restricted to the 
study of intermodal entanglement between Stokes mode 
and the vibration mode as they had considered a simpli- 
fied two-mode Hamiltonian. Thus the use of a completely 
quantum mechanical description of the Raman process, 
our solution, and the strategy to use more than one in- 
separability criterion have helped us to obtain a relatively 
more complete picture of the intermodal entanglement in 



the Raman processes. Further, if we look at the possi- 
bilities of different kinds of nonclassicalities summarized 
in Table U (see the rows corresponding to ac, be, and 
bd modes), then we can quickly recognize that the ex- 
istence of any one of the nonclassical phenomenon does 
not depend on the presence of the other. To be precise, 
entanglement, antibunching and squeezing are nonclassi- 
cal phenomena but they are independent of each other. 
However, Duan criterion in the present form implies in- 
termodal squeezing in one of the quadrature variable but 
the converse is not true. This fact can be clearly seen in 
the Table U where we note that except in bd mode inter- 
modal squeezing is possible in all other coupled modes. 
However, the Duan criterion of intermodal entanglement 
is satisfied only for the ac mode. 

In quantum optics, physical systems (matter-field in- 
teractions) are usually described by multi-mode bosonic 
Hamiltonians. The procedure followed in the present 
work may be applied directly to those systems to study 
the intermodal entanglement. It is expected that most of 
these systems will show intermodal entanglement. This 
is so because most of the quantum states are entangled. 
Separability is a very special case. A natural question 
should arise at this point: If entanglement is so common 
why are we looking for it? The answer lies in the fact 
that entanglement is one of the most important resources 
for quantum information processing and quantum com- 
munication but still it is not very easy to produce useful 
multi-partite entanglement. As many of the quantum op- 
tical systems described by bosonic Hamiltonian (includ- 
ing the system studied here) are experimentally achiev- 
able, useful intermodal entanglement may be produced 
by them. Entanglement obtained in such a system is 



expected to find application in controlled quantum tele- 
portation, quantum information splitting, dense coding, 
direct secured quantum communication etc. We conclude 
this paper with an optimistic view that the present work 
will motivate others to look for theoretical and experi- 
mental generation of useful multi-mode entanglement in 
other quantum optical systems. 
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Appendix A: Parameters for the solutions in Eq. (3) 



hi 
h 2 

h 3 

hs 
h e 
h 7 



exp(— iuj c t) 



„-iui c t 



vp — iu c t 



Aoj 2 

ye'*"":* 
Au> 2 



1] 
"I] 



i i(Au 1 +4 U2 )t_, 



AoJi J r Auj 2 



-1 



Au\ 



i/S.uj\t 



e A 

Aul 2 

,2 + „ — iuA 



ig 2 te 



= -h 5 



2 -iu> a t 
A c r iAuj 2 t 

Au,f [e 

2 -iu c t 

X c r iAuj 2 t 

Aa,2 L e 



1] 



l] 



Au 2 

ix 2 te -tu c t 



Auj 2 



fx 
h 

h 

h 

h 
h 
h 
h 



fji 

9z 

95 
96 



cxp(-iw a t), 

ge -iUat 



AU)! 



-iAuJit 



-i , 



XP r e iAcj 2 t 

Aoj 2 L 



1] 



Auji 

AW2 



i(Awi-/Vai 2 )t_l 



Awi — Aa;2 

i(A» I -Au 2 )t_j 



AaJ2 
„-*Aw,t 



Alui—Auj2 Awi 



Aw 2 



1] 



ig z te 



= h 



2 -iui a t 



Aw 2 



1] 



zx 2 ie~^* 



Acj 2 



= -/7 



exp(— 

-iu) b t 

^giAtJii 



Xge" 



AtJ2 (Awi — ALJ2) 



1], 

e i(Awi-Acj 2 )t _ ]1 



Au 2 (Awi +Aw2 ) 



2 —iui b t 

Aw 2 L 



rgiCAwx+Awa)* _ 

ig 2 te~ iMbt 



1] 



"55- 



(Al) 



(A2) 



Zi 

/2 



«5 
k 



exp(-iu)dt) 



-zAcJ2^ 



Aw 2 



1] 



Xge "" d ' Li(Aui-Au 2 )f _ ll 
AaJi(Awi-Acj 2 ) L J 
-iAuj 2 t 



h = 



Au)2 Acji 



1] 

3 --tAw 2 * _ ]1 



Aw 2 Acji 

Xge~'"d* r _{(A«i+A« 2 )i _ ll 
Awi(Acji+Aoj 2 ) I J 

xge~' 



Aw 2 



+ 



X e 



= Is 



9 



[1] O. Giihne, G. Toth, Phys. Rep. 474 (2009) 1. 

[2] A Peres, Phys. Rev. Lett. 77 (1996) 1413. 

[3] L M Duan, G Giedke, J I Cirac and P Zollar, Phys. Rev. 

Lett. 84 (2000) 2722. 
[4] H Hung and G S Agarwal, Phys. Rev. A, 49 (1994) 52. 
[5] J Lee, M S Kim and H Jeong, Phys. Rev. A, 62 (2000) 

032305. 

[6] R Simon, Phys. Rev. Lett. 84 (2000) 2726. 
[7] M Hillery and M S Zubairy, Phys. Rev. Lett. 96 (2006) 
050503. 

[8] M Hillery and M S Zubairy, Phys. Rev. A 74 (2006) 
032333. 

[9] M Hillery, H T Dung and H Zheng, Phys. Rev. A 81 
(2010) 062322. 

[10] G S Agarwal and A Biswas, New J. Phys. 7 (2005) 211. 
[11] C. H. Raymond Ooi, Q. Sun, M. Suhail Zubairy and M. 

O. Scully, Phys. Rev. A 75, 013820 (2007). 
[12] C. H. Raymond Ooi, Phys. Rev. A 76, 013809 (2007); 

Eyob A. Sete, and C. H. Raymond Ooi, Phys. Rev. A 

85, 063819 (2012). 
[13] A Miranowicz et al, Phys. Rev. A 82, 013824 (2010). 
[14] P Szlachetka, S Kielich, J Pefina and V Pefinova, J. Phys 

A: Math. Gen 12 (1979) 1921. 
[15] P Szlachetka, S Kielich, J Perina and V Pefinova, Opt. 

Acta 27 (1980) 1609. 



[16] J Perina, Quantum Statistics of Linear and Nonlinear 
Optical Phenomena (Kluwer, Dordrecht, 1991). 

[17] B Sen and S Mandal, J. Mod. Opt. 52 (2005) 1789. 

[18] B Sen, S Mandal and J Perina, J. Phys. B:At. Mol. Opt. 
Phys. 40 (2007) 1417. 

[19] B Sen and S Mandal, J. Mod. Opt. 55 (2008) 1697. 

[20] B Sen, V Pefinova, J Pefina, A Luks, and J Kfepelka, J. 
Phys. B: At. Mol. Opt. Phys. 44 (2011)105503. 

[21] P Grangier, Nature 438 (2005) 749. 

[22] L M Duan, M Lukin, J I Cirac and P Zoller, Nature 414 
(2005) 413. 

[23] V Pefinova, Luks, and J Kfepelka, J. Phys. A 44 (2011) 
035303. 

[24] S V Kuznetsov, O V Man'ko and N V Tcherniega, J. Opt. 
B: Quantum Semiclass. Opt. 5 (2003) S503-S512. 

[25] A Pathak, J Kfepelka and J Pefina , quant- 
ph/1210.3779vl. 

[26] D F Walls, Z. Phys. 237 (1970) 224. 

[27] Same values of x, g and \cti\ are used in the entire paper 
(unless otherwise specified). For spontaneous and par- 
tially spontaneous Raman proceses these values of |ai| 
are used used for non-zero |ai|'s. 

[28] S Sivakumar, J. Phys. B: At. Mol. Opt. Phys. 42 (2009) 
095502. 



10 




0.02 0.04 0.06 0.08 0.10 -25 



d) e) f) 

Figure 2: (Color online) Intermodal entanglement in stimulated Raman process by using HZ-1 criterion: The dotted 
line, dash-dotted line and the smooth line are used for the phase angle of the input complex amplitude ot\ for = 0, 
7rand 7r/2 respectively, a) intermodal entanglement is observed between Stokes mode and vibration-phonon mode for 

4> = 7r/2, b) intermodal entanglement is observed between pump mode and anti-Stokes mode, c) intermodal 
entanglement is observed between Stokes mode and anti-Stokes mode for (f> = and 7r/2, d) signature of intermodal 
entanglement is not observed between vibration-phonon mode and anti-Stokes mode, e) intermodal entanglement is 

observed between pump mode and vibration-phonon mode. 
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Figure 3: (Color online) Intermodal entanglement in stimulated Raman process by using HZ-2 criterion: The dotted 
line, dash-dotted line and the smooth line are used for the phase angle of the input complex amplitude ot\ for = 0, 

7r/2 and 7r respectively, a) intermodal entanglement is not observed between Stokes mode and pump mode b) 
intermodal entanglement is observed between Stokes mode and vibration-phonon mode, c) signature of intermodal 
entanglement is not observed between pump mode and anti-Stokes mode, d) intermodal entanglement is observed 
between Stokes mode and anti-Stokes mode only for <fi = 7r/2, e) intermodal entanglement is not observed between 
vibration-phonon mode and anti-Stokes mode, f) intermodal entanglement (time dependent) is observed between 

pump mode and vibration-phonon mode for <fi = tt/2. 
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Figure 4: (Color online) Intermodal entanglement in stimulated Raman process by using Duan criterion: The dotted 
line, dash-dotted line and the smooth line are used for the phase angle of the input complex amplitude u\ for = 0, 
7r/2 and 7r respectively. Intermodal entanglement is only observed between pump mode and anti-Stokes mode 

(FigS(c)). 



